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beaming pattern of photons produced by inverse Compton scattering. For an 
isotropic distribution of soft photons upscattered by nonthermal electrons with a 



By including Klein-Nishina effects, we generalize previous calculations of the 



power-law density distribution 71(7) oc 7~ p , embedded in a plasma moving with 
relativistic bulk speed we show that the observed radiation intensity is propor- 
tional to V 3+p , where T> is the Doppler boosting factor. This agrees with previous 
computations performed in the Thomson limit, where the observed spectral index 
is a = (p— 1)/2 and the beaming pattern J) i+2a . Independent of T>, Klein-Nishina 
effects limit the location of the peak energy e pea fcm e c 2 of the observed spectral 
energy distribution such that e pea k < l/eo> where e is the energy of the seed pho- 
tons in units of m e c 2 . Assuming that the seed photons originate in the broad line 
region, we demonstrate that the GeV emission of blazars is significantly modified 
by Klein-Nishina effects, the spectrum being softer than that calculated in the 
Thomson limit. We further show that the change in spectral index of the inverse 
Compton emission across e pea k can exceed the value of 0.5 predicted by compu- 
tations performed in the Thomson limit. The model spectra agree with OSSE 
and COMPTEL limits on this break without invoking the effects of differential 
absorption at the edge of a gamma-ray photosphere. 
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1. Introduction 

Inverse Compton scattering, a process in which a nonthermal distribution of electrons 
boosts the energy of a 'seed' or 'target' photon by a large factor, is commonly thought to 
be responsible for the production of gamma-ray photons in blazars (Ulrich, Maraschi, & 
Urry 1997). For very high energy electrons, the scattered photon can carry off a substantial 
fraction of the energy of the incoming electron. There have been several suggestions concern- 
ing the source of the target photons: optical/UV photons from an accretion disk (Dermer, 
Schlickeiser, & Mastichiadis 1992; Bottcher, Mause & Schlickeiser 1997), optical photons 
from the broad line region and infrared dust photons (Sikora, Begelman, & Rees 1994; 
Blazejowski et al. 2000). The radiation mechanism for such targets is usually called exter- 
nal Comptonization ('EC). Synchrotron photons produced in the jet can also act as seed 
photons for inverse Compton scattering, a process referred to as synchrotron self-Compton 
('SSC') scattering (Maraschi, Ghisellini, & Celotti 1992; Bloom & Marscher 1996; Inoue 
& Takahara 1996; Mastichiadis & Kirk 1997). The same synchrotron photons can serve as 
seed photons after being reflected back into the jet from the broad line clouds (Ghisellini & 
Madau 1996). In those blazars that have strong emission lines, the EGRET-detected GeV 
emission is probably dominated by inverse Compton scattering of the broad line photons 
(Sikora, Begelman, & Rees 1994; Sikora, et al. 1997). In this case, the plasma respon- 
sible for the GeV emission moves relativistically with a bulk Lorentz factor T through an 
approximately isotropic photon field. 

For such a 'blob' of plasma moving with velocity (3c = c(l — T -2 ) -1 / 2 , at an angle 9 
to the observer's line of sight, the beaming pattern of the observed power per unit solid 
angle and per unit frequency is T> 3+a for any process in which the emission is isotropic 
in the frame comoving with the blob (the 'blob frame'). Here, T> = l/[r(l — (3fi)} is the 
familiar Doppler factor with /i = cos 6, and a is the spectral index of the radiation. An 
isotropic electron distribution and a tangled magnetic field, produce, in the blob frame, 
isotropic synchrotron and SSC emission, which is, therefore, characterized by this beaming 
factor (Begelman, Blandford, & Rees 1984). The situation is different for inverse Compton 
scattering of photons on targets which are approximately isotropic not in the blob frame, 
but in the rest frame of the broad line region ('lab. frame'). In the Thomson limit, Dermer 
(1995) has shown that the beaming pattern in this case is X> 4+2a . In terms of the power- 
law index of the electron distribution, this is equivalent to £> 3+p . In § 2 we show that this 
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expression arises immediately from the transformation properties of the electron distribution 
and is valid for arbitrary electron energy. Thus, in the case of EC emission, the beaming 
is most simply expressed in terms of the underlying electron distribution, whereas for SSC 
and synchrotron emission it can be obtained straightforwardly from the observed photon 
distribution. In addition to the beaming pattern for a power-law distribution, Dermer (1995) 
also demonstrated that a characteristic energy, such as a cut-off or break in the spectrum, 
is observed at a photon energy that scales as e max oc T> 2 . This contrasts with the situation 
for synchrotron and SSC radiation, where the scaling is linear with T>. Once again, the 
difference can be understood from the transformation properties of electrons and photons, 
respectively. We show that the scaling changes when Klein-Nishina ('KN') effects intervene: 
for hard electron distributions (p < 3) with an upper cut-off, the peak in the spectral energy 
distribution is limited to a value e pea k < l/e that is independent of T>, where e is the energy 
of the seed photons in units of m e c 2 . (Energy units of m e c 2 are used throughout this work.) 

Assuming that the EGRET-detected GeV emission from blazars is due to inverse Comp- 
ton scattering of broad line region photons we show in § 3 that the spectra are significantly 
affected by KN effects, which result in a softening of the EGRET spectrum. This effect has 
been noted in connection with blazar spectra by Bottcher, Mause & Schlickeiser (1997), 
who computed model spectra using anisotropically distributed accretion disk photons as 
targets, rather than photons from the broad line region. We map the observed EGRET 
spectral indices to the index of the electron distribution, and find that these are grouped 
around the value p = 3.2, in agreement with recent theoretical work on particle accelera- 
tion in relativistic shocks (Bednarz & Ostrowski 1998; Kirk et al. 2000; Achterberg et al. 
2001). Synchrotron cooling or, alternatively, cooling by IC scattering in the Thomson limit 
produces a break of Ap = 1 in the electron spectrum. The spectrum of IC radiation in the 
Thomson limit therefore shows a break of Act = 0.5. However, some 7-ray sources that peak 
at MeV energies show a more pronounced break (McNaron-Brown et al. 1995), which has 
been interpreted as evidence for gamma-ray absorption by pair production (Blandford & 
Levinson 1995; Marcowith, Henri & Pelletier 1995). We show that spectral softening can 
explain these observations, provided the electron distribution is determined by synchrotron 
cooling. A summary of our results and their implications is presented in § 4. 



2. The beaming pattern of inverse Compton radiation 

Consider a blob of plasma moving relativistically with a bulk Lorentz factor T and 
velocity j3c, at an angle 9 to the observer's line of sight. In the frame of the blob the 
electrons are characterized by an isotropic power-law density distribution n' ('-/'), 
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n / (7 / ) = ^7 , - p J P(7i,72,7 / ), (1) 

Air 

where 7' is the Lorentz factor of the electron, A; is a constant, and ^(71,72,7) = 1 for 

7i < 7 < 72, and zero otherwise. Under the assumption that 7' 3> T, one can treat the 

electrons as a photon gas and make use of the Lorentz invariant quantity n/7 2 . The Lorentz 

factor 7 of an electron in the lab. frame is then 7 = 2>/ and the electron density 71(7) in 

the lab. frame is 

n( 7 ,/i) = n\i) (l\ = 1-^7-^(71^72^7)- (2) 

Given that the effective volume V e g of the blob is V e R = VT> (see Appendix A) , where V is 
the volume of the blob in the blob frame, the energy distribution of the effective number of 
electrons iV e flf(7,/i) is: 

kV 

iV cff ( 7 ,/i) = n( 7 ,/i)V; ff = — V 3+p 1 ' p P( ll V, l2 V, 1 ). (3) 

Consider now that this electron distribution will Compton-scatter seed photons of an 
arbitrary angular distribution. Since 7 > 1, we can make the approximation usually made 
in synchrotron theory that the outgoing photons move in the same direction as the scat- 
tering electron. Therefore, when observing under a certain angle only electrons moving in 
this direction contribute to the Compton luminosity. Since the effective number of these 
electrons is proportional to X> 3+p , the Compton specific luminosity (observed luminosity per 
energy interval per solid angle) is also proportional to D 3+p . Different seed photon angu- 
lar distributions will introduce an angle-dependent multiplication term in the calculation 
of the external Compton luminosity. An extreme case of a photon angular distribution is 
a monodirectional photon beam propagating along the direction of motion of the blob, for 
which Dermer, Schlickeiser, & Mastichiadis (1992) calculated the beaming factor of Comp- 
ton scattering to be D 3+P (l — /j,)( p+1 '' 2 . Note here that for 9 = 0° the inverse Compton 
luminosity is zero. It can easily be shown that for a monodirectional photon beam propagat- 
ing in the direction opposite to the direction of motion of the blob, as may be the case for 
the mirror model of Ghisellini & Madau (1996), the beaming factor is £> 3+p (l + jiY p+l ^ 2 . 

We now consider that the plasma blob is propagating through an environment permeated 
by an isotropic monoenergetic photon field of energy density U. The lab. frame rate of 
scatterings per final photon energy interval for an electron of Lorentz factor 7 is: 

Mp-*?llf( x ) (4) 

dtde ~ 4e o7 2/W ' [ ' 
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where <jt is the Thomson cross section. Jones (1968) introduced the 'head-on' approximation 
in which the seed photons are treated as coming from the direction opposite to the electron 
velocity. Using this, which is valid for for 7 ^> 1, and the full KN cross-section for inelastic 
Compton scattering, he showed that 

(4e o7 :r) 2 



m 



2x In x + x + 1 — 2x + 



2(1 
e 



4eo7a;) 



P(l/4 7 2 ,l,x), 



(5) 



x 



4e 7 2 (l - -) 



The maximum observed energy is 



^max,KN 



7 



46 072 2 ^ 2 



(6) 



1 + 4e 72^) ' 

In the case of Thomson scattering (7eo <C 1), Rybicki & Lightman (1979), assuming isotropic 
scattering in the electron frame, showed that 

2 



f(x) = -(l-x)P(l/AYA,x), 



x 



and that the maximum observed final energy is e ma x,T = 4eo7 2 £' 2 . 



(7) 



We now make the approximation that the outgoing photons are directed along the direc- 
tion of the scattering electrons, which is justified provided the electron angular distribution 
varies slowly over angular scales < I/7. To obtain the specific luminosity one integrates the 
scattering rate (4) over the electron energy distribution (3), and multiplies the result by the 
observed photon energy em e c 2 and by the photon number density n p = U/eom e c 2 



dL = p3+p 3WardJ_ e_ r°° ^ )f{x)p{ v v )rf (g) 

dealt lo7reo eo J\ 

In the Thomson case, for energies e m i nj T < e < e max ,T, where e m j ni T = 4eo7 2 £ >2 , the lower 
limit of the integration in equation (8) is , y m in = ( e /4eo) 1 ^ 2 , and the upper limit is 7 max = 72 T>. 
Performing the elementary integral using equation (7) we obtain: 



dL 

dedVt 



4e 



V 3 +P kVa T cU e_ ^ 
87re e 



(9) 



4e (3 + p)(72^) 2 l+p 

-(l+P)/2 o 



;i+p)(3+p) 
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For p > — 1 and e <C e max ,T we have r y m i n <C 7™^- Since the integrand is then steeper than 
7" 1 , the above result simplifies to 



dedVt 7re (l +p)(3 + p) \e 

The beaming of the observed radiation is the direct outcome of the electron beaming, and 
it is characterized by the electron index p. In the Thomson limit, the resulting spectrum 
is a simple power law with a spectral index a — (p — l)/2 and one can substitute for p in 
equation (10) to recover the X> 4+2a beaming pattern (equation 7 of D95). The spurious term 
(1 + / u) (a+1 ) in the result of D95, which, however, varies only slowly with jj, for viewing angles 
of interest, was introduced by the approximation that the seed photons in the frame of the 
blob are coming from a direction opposite to the direction of the velocity of the blob. In 
both D95 and here the maximum observed energy e maXt T, as well as any other characteristic 
energy, scale as oc V 2 , whereas in synchrotron and SSC they scale as T>. If, instead of 
observing at a fixed energy, we are interested in the specific luminosity measured at a break 
or cut-off in the spectrum, then the T> 2 scaling of the break energy introduces an additional 
V~ 2a factor, so that the specific luminosity at the break scales as T> 4 . The luminosity per 
logarithmic energy interval of the spectral feature, given by edL/dedQ, then scales as £> 6 , 
since e <xT> 2 . 

In the KN case, for energies e miny KN < e < e max ,KN, where e minjK N = 4e 7?£ )2 /( 1 + 
4e 7i£>), the lower limit of integration in equation (8) is found by setting x — 1 



In this case the integrand is also steeper that 7 -1 , and for 7 mi „ <^ 7 2 P =^ e <^ e max ,KN, 
the integration is dominated by the lower limit 7 m j„ which is independent of T>. Therefore, 
the beaming pattern "D 3+p is also valid in the general case of KN scattering. The maximum 
energy is given by equation (6), which is reduced to e maXt KN = 72^ when the high energy 
tail of the electron energy distribution is well into the KN regime, 72^>eo 3> 1, a behavior 
similar to that of synchrotron and SSC emission maximum observed energies. For electron 
indices p < 3, e pea k,KN cannot exceed significantly the energy at which KN effects become 
important and the scattering cannot be considered elastic. For a given seed photon energy eo 
this sets in for electrons with energies 7 ~ l/e^D. Setting this limiting value of 7 in equation 
(6) we obtain e pea k,KN ^ V e o> independent of V and 72, provided the system is well into the 
KN regime, 72^e ^> 1. 

We demonstrate these points in figure 1, where we plot the inverse Compton spectral 
energy distribution for three different observing angles for both the Thomson and KN cases 
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and for two different values of 72. We also plot the spectral energy distribution calculated 
using the approximation [a — &t f° r 7^ e o < 3/4, otherwise] (e.g. Chiaberge and Ghisellini 
1999, Blazejowski et al. 2000). In the Thomson case we use the analytical expression 
(10), while in the KN case we perform the integration in equation (8) numerically. The 
Thomson and the Klein Nishina distributions deviate from each other with the KN spectrum 
being softer. Note that the deviation is already significant at e m 10 4 , which corresponds 
approximately to electrons with Lorentz factor 7D m (e/eo) 1 / 2 ~ 4 x 10 4 in the lab. frame. 
Therefore, already at 'y'Deo ~ 0.2, the Thomson description is inadequate, and the KN 
formalism must be used. Both the maximum and peak energy of the Thomson spectral 
energy distribution scale as (72D) 2 . Contrary to this behavior, in the KN case the maximum 
energy scales as 72^, whereas the peak energy is insensitive to variations of both T> and 72 
and it is located at an energy e pea k,KN < l/ e o- The exact value of e pea k,KN is a function of 
the electron index p, with steeper electron power laws being characterized by lower e pea k^KN 
values. An increase in the upper cut-off 72 of the electron distribution by a factor of 10 affects 
only the steep high energy tail of the observed KN spectral energy distribution, leaving the 
peak energy and the peak luminosity unchanged. In general, as long as the scattering is 
KN limited ( 72^e 3> 1) , the peak energy will be insensitive to variations of both 72 and 
V, in contrast to the Thomson calculation and the synchrotron and SSC cases. The result 
based on the step function cross section [a = ot for 7^60 < 3/4, otherwise] is practically 
identical to the Thomson one, up to the cutoff energy l/eo- In the case of external Compton 
scattering of optical-UV photons in blazars, where l/e ~ 100 GeV, the spectrum calculated 
under this approximation in the EGRET regime is for practical purposes the same as that 
calculated in the Thomson regime. 

As shown in figure 1, the KN spectral energy distribution resulting from a power law 
electron energy distribution is not a power law, and one cannot assign a unique spectral 
index to it. The beaming pattern at a given energy is expressed through the electron index 
p and is not a simple function of the local spectral index — different parts of the spectrum 
have the same beaming pattern T> z+P independent of the local spectral index. This is in 
contrast to synchrotron and SSC emission from a power law electron distribution, where the 
beaming pattern depends only on the local spectral index, the steeper parts of the spectrum 
having a more pronounced beaming behavior. 



3. Applications to blazars 

The spectra of the EGRET -detected blazars (Hartman, et al. 1999) are described by 
simple power laws over the energy range 30 MeV - 10 GeV with no indication of a cut-off 



at high energy. The photon indices of those EGRET-detected blazars that display strong 
emission lines cluster around ~ 2.2 (Mukherjee, et al. 1997), indicating that the peak energy 
of the 7-ray spectral energy distribution in general lies at energies below the EGRET range. 
X-ray (e.g., Kubo et al. 1998) OSSE, and COMPTEL observations (McNaron-Brown et al. 
1995) confine this peak to between about 1 and 100 MeV. 

If the observed radiation is inverse Compton emission from optical/UV broad line seed 
photons (e ~ 10~ 5 ), a peak at ~ 10 MeV arises from electrons with 7Z? rs 1.5 x 10 3 . Since 
e^T* w 1.5 x 10~ 2 , the scattering can be adequately approximated by elastic Thomson 
scattering. Therefore, the peak at ~ 10 MeV is not connected to KN effects and must result 
from a break in the electron energy distribution. On the other hand, the EGRET-observed 
flux is KN affected, and cannot be described by elastic Thomson scattering. The 2 GeV flux 
(e?s4x 10 3 ) results from electrons with 7P rs 2 x 10 4 , corresponding to eo'-fV rs 0.2, a regime 
in which the KN steepening of the spectrum relative to the Thomson case is significant. 

In models in which particle acceleration competes with radiative losses and particle 
escape from the system, the electron energy distribution is characterized by 70, the Lorentz 
factor at which electrons are injected, 7&, the electron Lorentz factor at which the radiative 
cooling time equals to the escape time, and 7 max , the electron Lorentz factor at which 
the acceleration time equals the radiative cooling time. Between 70 and 7^ the electron 
distribution is a power law with index p, while above 75 the index steepens to p + 1, in 
the case where synchrotron cooling dominates. When the electron index p < 2, the peak 
of the inverse Compton spectral energy distribution is due to electrons with 7 rs '"fmax in 
the Thomson case, and the spectrum after the peak is expected to decrease abruptly. On 
the other hand, when the electron index p > 2, the peak of the inverse Compton spectral 
energy distribution is due to electrons with 7 rs 7^ in the Thomson case, and the spectrum 
after the peak is expected to follow a power-law behavior, up to a cut-off energy associated 
with '-/max. In the case of the EGRET-detected blazars the fact that after the Thomson 
dominated peak at rs 10 MeV there is a power law extension of the emission at least up 
to rs 10 GeV indicates that the peak of the spectral energy distribution is associated with 
electrons at 7 = 7^ where the radiative cooling time equals the escape time from the system, 
and that the electron index p > 2. 

If the EGRET-observed flux were due to Thomson scattering, the observed photon index 
s = a + 1 of a blazar would relate to the index p of the electron power law via p — 2s — 1. 
However, since KN effects steepen the EGRET spectrum, application of this relation results 
in an electron index that is too steep, as can be seen in the lower left panel of figure 2, where 
we plot the KN corrected photon index (solid line) as well as the relation s = (p + l)/2 
(dashed line). In this figure, the seed photon energy is eo = 5 x 10~ 5 , and the two-point 
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spectral index is calculated here and throughout the paper using the luminosity of the model 
spectrum at 100 MeV and 10 GeV. The maximum electron Lorentz factor is assumed to be 
72 22 > 1/eo) since the EGRET observations show no evidence of a cut-off. Under this 
assumption the observed photon index is not a function of 72 2X For photon indices s around 
2.2 the electron index that would result from Thomson scattering is steeper than the KN 
case by about 0.3. In the lower right panel of figure 2 we show a histogram of the observed 
photon indices of all blazars in the third EGRET catalog (Hartman, et al. 1999), except for 
sources with errors greater than 0.3 in their photon index determination and also excluding 
the three X-ray selected BL Lacertae (MKN 421, MKN 501, PKS 2155-304), which, because 
they lack line emission in their spectra, are believed to be predominantly SSC emitters, 

We see that the distribution peaks at s ~ 2.3 (a careful calculation for all the broad line 
emitting blazars gives s = 2.20 ± 0.05 (Mukherjee, et al. 1997)). For each of these sources 
we calculate the electron index p that would result in the observed photon index s, taking 
into account the correction due to the redshift of the source, and plot the results in the upper 
left panel of figure 2. The results cluster in a wide interval around p = 3.2 that includes the 
value p = 3.0. These values are expected from electron distributions with indices p = 2.2 
and p = 2.0, respectively, that steepen by Ap = 1 above 75. Since the EGRET range is 
at energies higher than the peak of the spectral energy distribution, the observed photon 
indices are compatible with acceleration at both relativistic (Kirk et al. 2000; Achterberg et 
al. 2001) and non-relativistic shocks, which predict p = 2.23 and p = 2 respectively for the 
uncooled electron distributions. Unfortunately, however, the scatter in the observed indices 
is of the same order as the difference between the predictions, so that the data do not allow 
us to distinguish between these two cases. 

Simple electron cooling considerations predict a spectral break of Aa = 0.5 in the 
transition before and after the peak energy e pea k of the observed energy distribution, as 
a result of the change in the electron index Ap = 1 for synchrotron dominated cooling. 
This appears to conflict with the combined OSSE, COMPTEL and EGRET measurements 
of some blazars, which found spectral breaks Aa > 0.5 (McNaron-Brown et al. 1995; 
Collmar et al. 1997). However, spectral breaks of Act > 0.5 between the COMPTEL 
and EGRET ranges are produced naturally for sources peaking at MeV energies, since the 
EGRET spectrum is softened by KN effects. We demonstrate this in figure 3, where we plot 
the spectral energy distribution due to inverse Compton scattering of optical seed photons 
by a broken power law electron distribution for both the KN (solid line) and Thomson 
calculation (broken line). The spectrum below the peak has a photon index s = (p+ l)/2 = 
1.6, since below the peak p = 2.2. Above the peak p = 3.2 and the Thomson spectrum has 
a a photon index s = 2.1, resulting in a break Aa = 0.5. The KN spectrum above the peak 
is steeper, and the two-point spectral index is calculated to be s = 2.30, which results in a 
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spectral break Aa = 0.70. 



4. Conclusions and discussion 

In this paper we generalize the calculation (D95) of the beaming pattern of inverse 
Compton radiation to the KN regime. In terms of the electron index, p, the beaming factor 
is T> 3+p , for arbitrary photon energy. The KN-limited spectrum produced by a power law 
electron distribution is curved, and the beaming pattern is not a simple function of the 
local spectral index as in the case of synchrotron emission or SSC. KN effects further limit 
the location of the peak of the spectral energy distribution to e pea k < l/ e o> where e is 
the seed photon energy. This peak energy is independent of the Doppler factor T> and the 
maximum Lorentz factor 72 of the electron distribution provided the scattering is KN limited 
(7 2 2}eo 3> 1), but does depend on the electron index, with e pea k shifting to lower energies 
as the electron distribution becomes steeper. If the GeV spectrum of the EGRET blazars 
with strong emission lines is due to external Compton scattering of optical seed photons 
(eo ~ 5 x 10~ 6 ), this component of the emission will steepen further due to KN effects at 
energies above 1/cq, corresponding to energies > 0.5 TeV. Thus, even in the absence of 
absorption due to the extragalactic IR background, the TeV flux of this EC component is 
negligibly small. 

In many cases, the spectra of blazars detected by EGRET and observed by OSSE 
and COMPTEL appear to peak in the region 1 - 100 MeV. Assuming that this emission 
is due to inverse Compton scattering of optical photons from the broad line region, the 
Thomson approximation is quite adequate at the peak energy. In this case, there should 
be a change in the spectral index measured on each side of the peak of Aa = 0.5, brought 
about by cooling via synchrotron or inverse Compton emission. In several cases, however, 
it appears that the spectra index determined from OSSE and COMPTEL measurements is 
harder than that in the EGRET-b&nd by Aa > 0.5. This apparent difficulty has led to 
the development of models in which the spectral break arises not from cooling, but from the 
differential absorption by pair production at the edge of a gamma-ray photosphere (Blandford 
& Levinson 1995; Marcowith, Henri & Pelletier 1995). 

Our results indicate, however, that the Thomson approximation is not adequate to de- 
scribe the spectrum in the EGRET range. Assuming that synchrotron radiation is the domi- 
nant cooling mechanism, we have shown that the spectral break between the OSSE / COMPTEL 
and EGRET ranges can substantially exceed 0.5, without invoking absorption due to pair 
production. Our computation assumes that the dominant cooling mechanism is synchrotron 
radiation. This is reasonable because the inverse Compton flux fjc is rarely observed to ex- 
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ceed the synchrotron flux f$ by more than a factor of 10. Application of the beaming formula 
with a — 1 implies for the ratio of the energy density Ue in target photons to the energy 
density U B in the magnetic field U E /U B = (fie/ fs)/V 2 . For V = 10 we have U E /U B < 0.1, 
so that synchrotron cooling dominates. We note that the issue of the large (Act > 0.5) spec- 
tral breaks in MeV blazars will be accessible to detailed observation following the launch of 
INTEGRAL. 

The spectral softening in the EGRET-b&nd brought about by KN effects implies un- 
derlying electron spectra which are ~ 0.3 harder than would be expected from a naive 
application of the formula p = (s — l)/2. As a result, the electron indices cluster around 
ps 3.2, which, allowing for cooling, corresponds to p Rs 2.2 below the break at 7^. This 
value is compatible with the predictions of shock acceleration theory at both at relativistic 
(p = 2.23) and nonrelativistic shocks (p = 2.0). However, the quality of the data do not 
permit us to distinguish between these two possibilities. 

The association of the peak of the spectral energy distribution with 7^ does not apply to 
all blazars. In X-ray bright objects, which include the TeV-emitting sources, the synchrotron 
spectral energy distribution is characterized by a sharp decline above the peak, as would 
be produced by a cut-off in the electron energy distribution, rather than a cooling break. 
Observations in the X-ray band of 'hard lags' (hard X-rays trailing in time behind soft X-rays) 
in MKN 421 (Maraschi, et al. 1999) suggest that at the peak energy the acceleration and loss 
time scales of the radiating particles are comparable (Kirk, Rieger, & Mastichiadis 1997). 
The corresponding electron spectrum is hard: p < 2. This difference in electron spectrum 
may reflect a fundamental difference in the underlying particle acceleration mechanism in 
these two distinct classes of objects. 

This work was supported by the European Union TMR programme under contract 
FMRX-CT98-0168. 



A. The effective observed volume 

Consider a blob containing a uniform distribution of electrons. In the rest frame of the 
blob, the differential number of electrons in interval d'j' moving within the solid angle Q' is 



dr/'dQ! Aiob 



Vn\i) (Al) 



-12- 



where V is the volume of the blob. In the lab. frame the blob moves along the x-axis at speed 
c/3, and the direction n to the observer makes an angle with this axis. The radiation flux 
F v observed is proportional to the product of the single particle emissivity and the electron 
distribution integrated over the emitting volume at the retarded time. This can be written 

p /"+00 

F v oc / d 3 x rftn(7,/i)5(t-x-n/c) (A2) 

./blob J — oo 

Apart from the argument of the delta function, the only dependence of the integrand on t 
arises from the limits of the integration over the moving blob, so that 

F v oc 71(7, //)Kflf = N et[ (j,/j,) (A3) 

where 



V P 



cff 



p p+00 

I d 3 x / dt5(t-x- h/c) (A4) 

./blob J— 00 



is the effective observed volume used in equation (3). To relate this to V, we rewrite Eq. (A4) 
in terms of the spatial coordinates in the blob frame, keeping the time t in the lab. frame. 
Since dx = dx'/T, dy' = dy and dz' = dz we have: 



hoc 



V cS = - d 3 x' I di<J (t - x • fi/c) 

r Jbiob 



which becomes, using x' = T(x — c(3t), 

V eS = ^ / d s x' I dtS (t - fiPt - fix'/(cT) - y/l - fi 2 y'/c) (A5) 

r Aiob j-00 ^ ' 

where we assume, without loss of generality, that n lies in the x-y plane. Performing the 
integral over t now yields 

K ff = _ l - I d 3 x' 



r(i - nP) Aiob 

VV (A6) 
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Fig. 1. — . The observed energy distribution due to inverse Compton scattering as a function 
of different observing angles for the KN (solid lines), the Thomson (broken lines), and the 
[ er = ot for 7X^0 < 3/4, otherwise] (dotted lines) cross sections for a blob of plasma that 
moves with a Lorentz factor T = 10 through an isotropic monoenergetic photon field. The 
seed photon energy is eo = 5 x 10~ 6 in units of m e c 2 , which corresponds to optical photons. 
The electrons in the blob frame are characterized by an isotropic power law distribution 
n( , ~f) oc 7 _p , p = 2.2, 71 < 7 < 72. For each case we plot in normalized units the result for 
both 72 = 10 5 and 72 = 10 6 , with the line corresponding to the higher 72 reaching higher 
photon energies. 
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Fig. 2. — Lower left panel: the electron spectral index p inferred from the inverse Compton 
scattering model as a function of the photon index s for the KN (solid line) and the Thomson 
(broken line) cases. Lower right panel: the distribution of the observed photon index s for 
blazars with strong emission lines taken from the third EGRET catalog (Hartman, et al. 
1999). Upper panel: the distribution of the electron indices p implied by the EGRET 
measurements, taking account of the individual source redshifts. 



-17- 








2 4 

log e 



6 8 



Fig. 3. — The observed energy distribution due to inverse Compton scattering for both the 
KN (solid line) and Thomson treatment (broken line) for a blob of plasma that moves with 
a Lorentz factor T = 10 through an isotropic monoenergetic photon field. The electrons in 
the blob frame are characterized by an isotropic broken power law distribution 71(7) oc 7~ p , 
with p = 2.2 for 10 < 7 < 2 x 10 2 , and p = 3.2 for 2 x 10 2 < 7 < 10 5 . We plot, in 
normalized units, the energy distribution due to inverse Compton scattering observed at an 
angle 6 = 1/T. The shaded area corresponds to the EGRET range of observation. 



